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ABSTRACT 

We determine an absolute calibration of stellar mass-to-light ratios for the densest 
~ 3% of early-type galaxies in the local universe (redshift z ~ 0.08) from SDSS 
DR7. This sample of ^ 4000 galaxies has, assuming a Chabrier IMF, effective stellar 
surface densities, So > 25OOM0/pc^, stellar population synthesis (SPS) stellar masses 
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logio(Mspg/M0) < 10.8, and aperture velocity dispersions of (Tap = 1681^4 km s-i 
(68% range). In contrast to typical early- type galaxies, we show that these dense early- 
type galaxies follow the virial fundamental plane, which suggests that mass-foUows- 
light. With the additional assumption that any dark matter does not follow the light, 
the dynamical masses of dense galaxies provide a direct measurement of stellar masses. 
Our dynamical masses (Mdyn), obtained from the spherical Jeans equations, are only 
weakly sensitive to the choice of anisotropy (/3) due to the relatively large aperture 
of the SDSS fiber for these galaxies: i?ap — l-5i?c- Assuming isotropic orbits {P = 0) 
we find a median log;iQ(Mdyn/Msps) = 0.233 ±0.003, consistent with a Salpeter IMF, 
while more bottom heavy IMFs and standard Milky- Way IMFs are strongly disfavored. 
Our results are consistent with, but do not require, a dependence of the IMF on 
dynamical mass or velocity dispersion. We find evidence for a color dependence to the 
IMF such that redder galaxies have heavier IMFs with Mdyn/Msps oc {g — 7.)i i3±0 09^ 
This may reflect a more fundamental dependence of the IMF on the age or metallicity 
of a stellar population, or the density at which the stars formed. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters 
- galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

The stellar initial mass function (IMF) is a fundamen- 
tal property of a stellar population, with wide ranging 
implications for many areas of astrophysics. Observations 
in the Galactic disk suggest that the IMF has a power- 
law shape, dN/dm oc , with x ~ —2.35 at masses 
above m ~ IMq (Salpeter 1955), and that it turns over 
at lower masses (Kroupa 2001; Chabrier 2003). For sim- 
plicity, the IMF is generally assumed to be universal, al- 
though there is increasing evidence that the IMF depends 
on the mass (or velocity) of a galaxy. Spiral galaxies require 
IMFs lighter than Salpeter (Bell & de Jong 2001; Button 
et al. 2011 a,b; Brewer et al. 2012) and maybe even lighter 
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than Chabrier (Bershady et al. 2011). While elliptical galax- 
ies favor Salpeter-type IMFs (Auger et al. 2010), and pos- 
sibly bottom heavy {x ~ —3) IMFs in the most massive 
ellipticals (van Dokkum & Conroy 2010). 

Total mass measurements from dynamics or strong lens- 
ing only give an upper limit to stellar mass-to-light ratios 
due to the unknown dark matter fraction. However, if we 
can identify a population of galaxies that are baryon domi- 
nated (i.e., the dark matter fraction within a specific radius 
< 10%), these can be used to provide a direct dynamical 
constraint to stellar mass-to-light ratios, and hence the IMF. 

In this Letter we show that early-type galaxies with 
the highest stellar surface densities within an efi'ective ra- 
dius have a Fundamental Plane (FP; Dressier et al. 1987; 
Djorgovski & Davis 1987) correlation that is consistent with 
the virial prediction: At fixed SPS stellar mass, velocity dis- 
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strong evidence that the dynamical masses are proportional 
to the stellar masses. Under the additional assumption that 
any dark matter does not trace the light, the dynamical 
masses are equal to the stellar masses. Here we use spheri- 
cal Jeans models to calculate dynamical masses and compare 
these with SPS stellar masses. We adopt a cosmology with 
Q.A = 0.7, fim = 0.3 and Hq = 70 kms"^ Mpc"\ 

2 SAMPLE SELECTION 

Our selection of early-type galaxies is designed to select 
non-star forming galaxies regardless of their morphology 
or structure, and is the same as in Button et al. (2Qlla). 
Our starting sample consists of ~ 685000 galaxies from the 
Sloan Digital Sky Survey (SDSS; York et al. 2000) data 
release seven (DR7) (Abazajian et al. 2009) with spectro- 
scopic redshifts (0.005 ^ z; ^ 0.2), stellar velocity dis- 
persions, (Jap, from the "Princeton" catalog, structural pa- 
rameters (using the Sersic n = A plus n = 1 fits) from 
Simard et al. (2011), and stellar masses computed assuming 
a Chabrier (2003) IMF from the MPA/JHU group (available 
at http://www.mpa-garching.mpg.de/SDSS/DR7/). We se- 
lect galaxies that have been spectroscopically classified as 
early-type (eCLASS < 0), have red {g — r) colors, and 
have minor-to-major axis ratios greater than 0.5. We also 
apply redshift dependent minimum stellar mass. Our se- 
lection yields ~ 140000 early-type galaxies at redshifts of 
z ~ 0.079t2 g^J (68% range). 

Using the stellar masses, Msps, and the circularized 
r— band effective radii, i?c, we compute the average sur- 
face density within the effective radius: Ec = Msps/(27ri?^)- 
We select a sample of dense galaxies with the criteria Ee > 
2500Mqpc2. This yields ~ 5000 galaxies, i.e., ~ 3.5% of our 
full early-type galaxy sample. While the exact value of this 
minimum density is somewhat arbitrary, its approximate 
value will become apparent below. In order to minimize the 
impact of galaxies scattered to high stellar densities by the 
~ 0.1 dex measurement errors on stellar masses, we clean 
the sample by using the relation between Msps/Lr and a"ap. 
This relation has a slope of ~ 0.4 and a scatter of ~ 0.06 
dex. We remove the outlier galaxies that are more than ±3(j 
offset from the median relation. This removes ~ 400 galaxies 
leaving a dense sample of ~ 4600. 

The velocity-mass and size-mass scaling relations for 
our full and dense samples of early-type galaxies are shown 
in Fig.[T] The median relations of the full samples are shown 
with points and solid lines, while the high density galaxies 
are shown with a color scale. In the size-mass relation our 
selection corresponds to a line of slope 1/2. All of the dense 
galaxies are offset by at least la from the median size-mass 
relation. In the velocity-mass plane the dense galaxies have 
a slope of 0.25 and a are offset to high velocity dispersions. 

3 EVIDENCE FOR MASS FOLLOWS LIGHT 

The assumption of mass-foUows-light (MFL) is sometimes 
used when modeling the stellar kinematics of luminous early- 
type galaxies (e.g., Cappellari et al. 2006). There are two 
issues with MFL models: 1) MFL models only provide an 
upper limit to the stellar mass due to the unknown dark 
matter fraction; 2) If there is a significant amount of dark 
matter, then the assumption of MFL likely fails and the 
dynamical masses will be biased. One can get around these 
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Figure 1. Size vs stellar mass and velocity dispersion vs stellar 
mass relations for our full and high-density samples of early-type 
galaxies. The relations for the full sample are shown with lines: 
median (solid); 16th and 84th percentiles (dashed); 2.5th and 
97.5th percentiles (dotted). The color shading shows the relations 
for our high-density galaxies. In the size-mass relation the density 
selection is given by the red line with slope 1/2. 

two problems by selecting galaxies for which the dark matter 
fraction is negligible. To select such a sample we use the 
correlations between velocity and size at fixed stellar mass. 

The simplest mass follows light models have the scaling 
cr^(ii) oc Mstur/R- It is well known that the Fundamental 
Plane is tilted with respect to this model. Here we express 
the FP relation as the correlation between residuals from the 
velocity dispersion vs stellar mass (A log (Tap) and effective 
radius vs stellar mass (A log Rc) relations. We label the slope 
of the correlation as (?vr = A log (Tap / A log R^. The value of 
cVr is related to the dark matter fraction within the radius 
the velocities are measured (Courteau & Rix 1999; Dutton 
et al. 2007). Models with no dark matter nominally have 
(JvR = ^1/2, while dark matter dominated galaxies have 
9vR > 0. 

In Fig. [5] the solid black line shows the FP relation for 
all early-type galaxies. Fitting a linear relation gives a slope 
9vR ~ —0.300 ± 0.001 (red line) which is shallower than 
the virial prediction (dotted line). This recovers the well 
known tilt of the FP. It is also apparent that the local value 
of dvn varies monotonically with size offset. Galaxies with 
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Figure 2. Correlation between residuals of the Cap — A^sPS and 
Re — Afgps relations. The virial fundamental plane has a slope 
of —1/2 (dotted line). The observations are given by the black 
line. For i?e — M^SPS residuals within ±0.2 dex (which is 95% of 
galaxies) the relation is linear with A log (Tap = — 0.30 A log i?c, 
recovering the well known tilt of the fundamental plane. However, 
for larger size offsets the data clearly depart from a linear relation, 
representing curvature to the fundamental plane. 



the most negative offsets have i9vr — —0.5, while galaxies 
with the most positive offsets have dyn — 0.0. This trend 
is qualitatively consistent with the expectations for galaxies 
embedded in extended dark matter haloes. 

From Fig. [2] we see that galaxies with Alogi?e —0.2 
have a local slope of i9vr — —0.5 and thus are consistent with 
the assumption of MFL. This offset is the motivation behind 
our particular density selection as shown in Fig. [1] The col- 
ored points in Fig. [3] show the FP relation for galaxies in our 
dense sample, split into high and low stellar masses. The low 
mass galaxies in our sample, 9.6 < \og-^Q[Msta.i / Mq) ^ 10.8, 
(red points) have 9vr — —0.5 which suggests that these 
galajdes are dominated by stars within an effective radius 
(Courteau & Rix 1999). However, the massive galaxies, 
10.8 < logjo(Mstar/M0) < 11.6, (blue points) have a weaker 
correlation, suggesting that mass does not follow light. 

The possibility that dark matter follows the light in 
the inner regions of elliptical galaxies has been discussed by 
Thomas et al. (2011). However, this scenario seems unlikely 
for our dense galaxies, since at the effective radius the slopes 
of the stellar mass profiles are d In pstar/d Inr ~ —2.5, while 
the dark matter is expected to have dlnpDM/dlnr > —2. 

Thus in what follows we can only strictly justify MFL 
models for the lower mass dense galaxies. However, using re- 
lations from Maccio et al. (2008) and Button et al. (2010), 
we estimate the dark matter fractions within i?c to be just 
~ 5%. We will still compute MFL models for the high mass 
dense galaxies, with the caveat that dark matter may pro- 
vide a non- negligible contribution to the dynamical masses. 
Our dense-low mass sample consists of ~ 4000 galaxies with 
median SPS masses Msps ~ 3 x 10^° M©, velocity dis- 
persions (Tap = 168^34 km s~^, and fiber aperture radii of 
i?ap — 1.5i?c- The relatively large aperture radii results in 
an additional benefit when deriving dynamical masses from 
aperture velocity dispersions. For small apertures relative 
to the effective radius (e.g., one eighth of an effective ra- 
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Figure 3. Correlation between residuals of the (Tap — Afgpg and 
Rc — Mgps relations for dense galaxies. For reference, the solid 
line shows the data for all galaxies from Fig. [J] but zoomed 
into the upper left quadrant (A log (Tap > 0, Alogi?e < 0). 
The red and blue points show the relation for dense galaxies 
(So > 2500 Mqpc-^) split into high and low masses at Mgpg = 
IO^^'^Mq. The low mass galaxies (red points) have a slope equal 
to the virial fundamental plane which suggests that they are 
baryon dominated within an effective radius. 

dins), the derived dynamical masses are strongly dependent 
on the anisotropy profile. When the aperture is larger than 
the effective radius the anisotropy is only of secondary im- 
portance. 

4 DYNAMICAL MASS VS STELLAR 
POPULATION SYNTHESIS MASS 

For each galaxy in our high-density sample we calculate 
a dynamical mass using the following procedure. Given a 
galaxy light profile (parameterized by circularized and de- 
projected Sersic n = 4 and n — 1 components — see Button 
et al. 2011a) with mass normalized to Msps, we calculate 
the radial velocity dispersion profile by solving the spheri- 
cal Jeans equations assuming a constant anisotropy, /3 (e.g., 
Binney & Mamon 1982). We then compute the projected 
velocity dispersion (convolved with 1.4 arcsec FWHM see- 
ing) within the 3 arcsec diameter aperture used by SBSS. 
Finally, we compute a dynamical mass by scaling the model 
aperture velocity dispersion to match the observed aperture 
velocity dispersion: Mdyn = Msps(o-ap,obs/o-ap,modci)^- 

Note that since the model velocity dispersions scale as 
Mgp| (because we assume mass- follows- light) our dynami- 
cal masses are actually independent of the SPS masses. Our 
procedure also takes into account structural non-homology, 
which can significantly effect dynamical masses (e.g., Taylor 
et al. 2010). While our models do not account for rotation, 
Cappellari et al. (2006) have shown that two and three inte- 
gral models yield dynamical masses that are consistent with 
those obtained from a simple virial relation Mdyn cx: a^Re 
for both fast- and slow-rotators, which thus supports our use 
of spherical Jeans models. In addition we have verified that 
the dynamical masses show no correlation with axis ratio. 

For isotropic orbits (/? = 0) the median offset between 
dynamical and SPS masses for the full dense sample is 
AiMF = log(Mdyn/Msps) = 0.234 ± 0.003. For the low- and 
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Figure 4. Relation between dynamical mass, M^y^, assuming 
/3 = 0, and SPS stellar mass, Msps, assuming a Chabrier (2003) 
IMF, for dense galaxies. The relation between dynamical and SPS 
masses has a slope of 1.000 it 0.007 (for a forward fit, solid line), 
and 1.225 it 0.009 (for an inverse fit, dashed line). The zero point 
offset is in excellent agreement with a Salpeter IMF. The dynam- 
ical masses are, on average, inconsistent with lighter Milky- Way 
type IMFs (Chabrier), or bottom heavy {x = 3.0, x = 3.5) IMFs. 

high-mass subsamples we find Aimf = 0.233 ± 0.003, and 
AiMF = 0.236 ± 0.005, respectively. As expected the dy- 
namical masses are only weakly dependent on the assumed 
velocity anisotropy For reasonable values of velocity 
anisotropy (0.0 < /3 0.5) for early-type galaxies (Gerhard 
et al. 2001; Koopmans et al. 2009) the dynamical masses dif- 
fer by just 4% (with lower dynamical masses for higher /3). 
Thus, on average, our results are in excellent agreement with 
a Salpeter IMF (Aimf — 0.23), but are strongly inconsistent 
with lighter Milky- Way IMFs (e.g., Kroupa 2001; Chabrier 
2003), or heavier IMFs found in the centers of more massive 
ellipticals (van Dokkum & Conroy 2010). 

4.1 Variation of the IMF with galaxy mass 

Fig.|4]shows the relation between dynamical and SPS masses 
for /3 = 0. The dotted lines show the relations expected for 
various IMFs, ranging from Chabrier to a; = 3.5. Due to 
the non-uniform distribution of Msps (which is well approx- 
imated by a log-normal), and the scatter in Mdy-n/ Msps of 
~ 0.15 dex, the slope of the true relation between Mdyn and 
Msps will depend on the origin of the scatter. If the scat- 
ter is entirely due to Mdyn the slope is 1.000 ± 0.007 (solid 
line), while if the scatter is entirely due to Afsps the slope 
is 1.225 ± 0.009 (dashed line). In practice there are likely 
errors on both masses, and the true slope of the relation 
between Mdyn and Msps will fall between these two limit- 
ing cases. The nominal measurement errors on Msps of 0.09 
dex, imply a true relation of Mdyn oc Mgps . 

Since dynamical mass is proportional to the square of 
velocity dispersion, our results are consistent with the no- 
tion that the IMF depends on the velocity dispersion of 
the galaxy, with heavier IMFs in higher dispersion galax- 
ies (Treu et al. 2010; van Dokkum & Conroy 2011). Such 
a trend is also consistent with the lighter IMFs favored for 
spiral galaxies (e.g., Bershady et al. 2011), which tend to 
have lower velocities than ellipticals. However, since veloc- 
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Figure 5. Dependence of the IMF offset parameter, Aimf = 
l°gio(^dyn/Msps)i On (g — r) galaxy color for low-mass dense 
early- types. Redder galaxies have heavier IMFs. 

ity dispersion correlates with other physical galaxy proper- 
ties (e.g.. Graves et al. 2009), variation in the IMF may be 
driven by other properties of a stellar population, such as 
age or metallicity. 

4.2 Variation of the IMF with galaxy color 

Fig. [5] shows that the IMF offset parameter is correlated 
with color, with the sense that redder galaxies have heavier 
IMFs. This relation is given by 

Aimf = 0.222(±0.004) + 1.12(±0.09)[(g - r) - 0.75]. (1) 

Here the colors and luminosities have been k-corrected to 
z = (but not evolution corrected). Extrapolating Eq.[l]to 
bluer colors, a Chabrier IMF (i.e., Aimf = 0) occurs at a 
color of ((? — r) = 0.56, which is typical for a Msps ~ 3 x 
W^°Mq late-type galaxy (e.g., Fig.l of Dutton et al. 2011a). 

A color dependence to the IMF could plausibly be re- 
lated to the idea that the characteristic mass scale of the 
IMF is set by the Jeans mass of the ISM, which depends on 
both temperature and density (Bate & Bonnel 2005; Larson 
2005) . A density dependence would take the form of heavier 
IMFs for stars formed in denser regions. Since more massive 
galaxies formed their stars at higher densities, and more 
massive galaxies have older and redder stellar populations, 
it follows that galaxies that formed their stars at higher den- 
sities (with heavier IMFs) should have redder colors. 

4.3 Dynamical mass-to-light ratio vs color 

With the goal of providing a useful estimator of stellar mass- 
to-light ratios, Fig.[S]shows the relations between dynamical 
(and SPS) mass-to-light ratio vs {g — r) color for a model 
with isotropic orbits (/3 = 0). SPS mass-to- light ratios scale 
as Msps/Lg oc (fl-r)^'®^ and Msps/ir oc (5-r)^'^\ These 
slopes are very similar to those of Bell et al. (2003), who 
report 1.52 and 1.10 for the g- and r-bands, respectively. 
The relations between dynamical mass-to-light ratios and 
color in the g- and r-bands from Fig. [6] are given by the 
following equations: 

Mdyn/is = 0.640(±0.003) + 2.61(±0.09)[(.g - r) - 0.75], (2) 
Mdyn/ir = 0.532(±0.003) + 2.21(±0.09)[(3 - r) - 0.75]. (3) 
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Figure 6. Relations between mass-to-light ratios and (g — r) 
color for a model with isotropic orbits (/3 = 0). Dynamical mass- 
to-light ratios (Mjyn/L) have stronger scalings with color than 
SPS mass-to- light ratios (Mspg/L). 

The steeper slopes of the M^yn/L vs. (g — r) compared to the 
Msps /L vs {g — r) relations confirm that the trend between 
IMF offset and galaxy color shown in Fig. [5] is not due to 
correlations between Msps and (g — r). 

5 SUMMARY 

We have determined an absolute calibration of stellar mass- 
to-light ratios from the densest ~ 3% of early-type galax- 
ies in the local universe (redshifts z ~ 0.08) from the 
SDSS DR7. The fundamental plane scaling relations of these 
galaxies (with SPS masses logjQ(Afsps/M0) < 10.8, assum- 
ing a Chabrier IMF) are consistent with the virial plane: 
(Jap oc Msps/Rc, suggesting that mass-foUows-light. With 
the additional assumption that any dark matter will not 
follow the light, we infer that these galaxies are dominated 
by stars, i.e., they contain little or no dark matter within an 
effective radius. 

We compute dynamical masses using spherical Jeans 
models with constant velocity anisotropy, /3. For isotropic 
orbits (/3 = 0) the dynamical masses are on average 0.233 ± 
0.003 dex heavier than the SPS masses obtained with a 
Chabrier (2003) IMF. For radially anisotropic orbits with 
P = 0.5 the median offset is just 0.02 dex lower. The depen- 
dence on the anisotropy is small because the velocity disper- 
sions are measured within a fairly large aperture: ~ l.Siic. 

The relation between dynamical and SPS masses, 
M^yn oc M^pg, has a slope 1.00 ^ a ^ 1.23 depending 
on the source of the scatter in this relation, with a most 
likely value of a = 1.1. Our results are thus consistent with, 
but do not require, a dependence of the IMF on dynamical 
mass or velocity dispersion. We find a strong dependence 
between Mdyn/Msps and {g — r) color, with the sense that 
redder galaxies have heavier IMFs. If this trend refiects a 
more fundamental dependence of the IMF on the age and/or 
metallicity of a stellar population, it implies that the IMF 
varies within as well as between galaxies. 

Future applications of this method. This method has the 
potential to be applied to galaxies at higher redshifts, and 
thus probe the evolution of the IMF. Since the sizes of early- 
type galaxies are known to decrease as we look backward in 



time (e.g., Trujillo et al. 2006), we expect our method to 
be applicable to a larger fraction of early-type galaxies, and 
specifically higher mass galaxies, at higher redshifts. 
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